Higher-order mesoscopic fluctuations in quantum wires: 
Conductance and current cumulants 
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We study conductance cumulants {((?")) and current cumulants Cj related to heat and electrical 
transport in coherent mesoscopic quantum wires near the diffusive regime. We consider the asymp- 
totic behavior in the limit where the number of channels and the length of the wire in the units of 
the mean free path are large but the bare conductance is fixed. A recursion equation unifying the 
descriptions of the standard and Bogoliubov-de Gennes (BdG) symmetry classes is presented. We 
give values and come up with a novel scaling form for the higher-order conductance cumulants. In 
the BdG wires, in the presence of time-reversal symmetry, for the cumulants higher than the second 
it is found that there may be only contributions which depend nonanalytically on the wire length. 
This indicates that diagrammatic or semiclassical pictures do not adequately describe higher-order 
spectral correlations. Moreover, we obtain the weak-localization corrections to Cj with j < 10. 

PACS numbers; 73.63.Nm, 74.25.Fy, 72.15.Rn, 05.40.-a 



I. INTRODUCTION 

In a quantum wire (quasi-one-dimensional geometry), 
the quantities of interest, such as conductance, are not 
selfaveraging.^ Thus, a statistical description of a vari- 
ety of quantum effects, such as universal conductance 
fluctuations, requires the conductance cumulants ((g")) 
or conductance distribution P{g) to be studied. Cumu- 
lants of current Cj , like C2 that is directly related to shot 
noise, also yield information which is not contained in the 
averaged conductance. Higher-order mesoscopic fluctua- 
tions, though smaller, provide important additional in- 
formation about the transport process. For example, the 
recent experiments^ where an anomalously asymmetric 
P{g) was detected have spurred dispute"^ and interest 
on higher conductance cumulants in mesoscopic wires. 
This is due to the fact that noninteracting one-parameter 
scaling models suggest that near the diffusive region the 
conductance distribution P{g) is close, but not identical, 
to Gaussian shape.'' For the moment, current-cumulant 
measurements on the third cumulant C3 have been car- 
ried out^ and a scheme for a detection of C4 has recently 
been put forward.^ 

The fundamental symmetries of the Hamiltonian of a 
disordered conductor manifest themselves in the statisti- 
cal properties of the energy levels and particle states. 
The implications of symmetry are conveniently ana- 
lyzed within the framework of random matrix theoryii 
There are altogether ten different random matrix theo- 
ries for different symmetric spaces^ of the Hamiltonian.^ 
The symmetry classes are customarily referred by the 
Cartan's symbol for the symmetric space of the corre- 
sponding Hamiltonian.^ In this paper we consider seven 
symmetry classes: A(I,II), C(I), and D(III). Disordered 
normal metals may be studied through standard or 
Wigner-Dyson (WD) universality classes A(I,II)iiS The 
BdG classes C(I) and D(III) are appropriate, e.g., for 
"disorder-facilitated" quasiparticle transportii in uncon- 
ventional superconductors. Compared to the WD en- 



sembles, an extra degree of freedom arises for the BdG 
classes, since a distinction between particlelikc and hole- 
like quasiparticles of the BdG formalism is made. Fur- 
ther, the systems are classified according to the pres- 
ence or absence of time-reversal (TR) and spin-rotation 
(SR) symmetry. The corresponding symmetric spaces 
are characterized by the multiplicities of the ordinary 
and long rootsi^ denoted toq and nii (see, e.g.. Table 
1 of Ref. IT^ . In addition to the two symmetry pa- 
rameters, the degeneracy d of transmission and reflec- 
tion eigenvalues has to be taken into account. Besides 
having definite values for certain symmetries, mg and 
mi may also be considered as interpolation parameters 
between different symmetry classes. In the absence of 
electron-electron interactions, phase randomization and 
local maximum entropy principles imply in a quasi-one- 
dimensional geometry a scaling equation, the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equationii The DMPK 
equation yields the evolution of P{g) as a function of di- 
mensionless length s. We set s — L/N£, where L is the 
length of the wire, N is the number of channels, and i 
is the elastic mean free path. In a wire geometry, the 
DMPK equation has also been generalized for the BdG 
structuresi^ and for the so-called chiral classes 

Near metallic regime the most important mechanisms 
inducing deviations from Gaussian P{g) are weak local- 
ization (WL) and weak antilocalization (WAL). WL and 
WAL result from the interference of the closed trajecto- 
ries of the electron and their time-reversed conjugates. 
Constructive interference (WL) leads to a decrease of 
g while destructive interference (WAL) enhances con- 
ductance. In a conductor with length shorter or com- 
parable with the coherence length, in the absence of 
magnetic field, destructive interference may be realized 
in a material with strong spin-orbit scattering whereas 
spin-rotation symmetry leads to constructive interfer- 
ence. The W(A)L corrections to conductivity have been 
used to study experimentally the quantum coherence in 
mesoscopic structures (see, e.g., Ref. ITtI for recent mea- 



surements). It is very difficult to evaluate the exact influ- 
ence of W(A)L on the values of (((?")) with n > 3 by such 
field theoretic approaches as diagrammatic techniques^ii^ 
or by the nonlinear sigma model^ 

A nonperturbative treatment of the first and second 
cumulants based on Fourier analysis on the supersym- 
metric manifold was given in Ref.Mfor WD classes. The 
behavior of the lowest conductance cumulants ((.g^'^)) has 
been intensively computed for the WD ensembles (to/ 
1) recently, especially in the metal-insulator crossover re- 
gion (see, e.g., Refs. I21I andl23). In Ref. 0, based on a 
2 + e expansion, the expression 



II. CUMULANTS 



A. Method 



The starting point for our analysis is the generalized 
DMPK equation which leadsi^^ 
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was presented for the standard unitary ensemble 
(too = 1, To; = 1). Likewise, for the symmetry class AI, 
but for a quasi-one-dimensional conductor (where one has 
(g) fa 1/s), the same equation was put forward in Ref.l23l 
This formula has been widely accepted until now. 

For the BdG wires near metallic region, the depen- 
dence of {g) on Too and to/ may be found in Ref. ITsl 
Leading terms of ((5^)) (universal conductance fluctua- 
tions) and some correction terms for ((g^)) and for C2 
may be found in Ref. |2j. An essential singularity has 
been shown to occur for the first and second cumulants 
in the small s expansion in the chiral unitary classs^^^ 
Macedo considered in Ref. i2Q the cases with too = 2 and 
< mi < 2 and found that the third cumulant contains 
only a component which is nonanalytic in s at s = 0. Lo- 
calization in superconducting wires has been discussed in 
Refs. [3 and [23. The conductance cumulants with rt < 4 
were computed in Ref. '21' for the WD classes A and AI 
by a Monte Carlo method. 

In this paper we present a recursion equation unifying 
the description of the WD and BdG symmetry classes 
and yielding the cumulants of the order n < 1/s. For 
the two BdG classes with TR symmetry we find that 
the higher-order cumulants {n > 3) contain no contribu- 
tions that are analytic in s at s = 0. We elucidate the 
dependence of the higher cumulants on the universality 
class and give values for ((g")) in terms of too and to;. 
We emphasize that even though Eq. |^ is correct for 
n = 1,2, for 2 < n < 1/s there exists a more appropri- 
ate expression, our Eq. (|12|l . Furthermore, we calculate 
the weak-localization corrections to the current cumu- 
lants Cj with j < 10. We consider a disordered quan- 
tum wire, i.e., a quasi-one-dimensional geometry. For the 
BdG universality classes we study heat transport whereas 
for the WD classes our results apply also for electrical 
transport. We calculate the cumulants at zero temper- 
ature, zero frequency, and at low voltage in the limit 
N — > 00, L/£ 00, s = constant, near the diffusive re- 
gion, where one has 1/A^ ^ s ^ 1. 



The variables {Xi}^i are related to the transmission 
probabilities {Ti}fLi of the channels {i}fLi through 
Ti = (1 + Xi)~^ while Ws{X) is the distribution function 
for A = (Ai, . . . , Aat). The Jacobian JmoW is given by 
JmoW = I\i<j - Ajl™". In a long wire, for the sym- 
metry classes D(III), additional terms enter the DMPK 
equationiSl Near diffusive regime such components, how- 
ever, are irrelevant and we will ignore them. 

In order to calculate the cumulants we adopt a method 
reminiscent to the moment expansion method introduced 
m Refs. Inland (291 It is convenient to introduce the mo- 
ment generating function Zs (q) and the cumulant gener- 
ating fmiction (CGF) fs{<i), 



lnZs{q) = ln(cxp(-q • T))^ = (^s(q) 



(3) 



in order to systematically solve the DMPK equation. 
Here one has T = (Ti, . . . , Ta,), Ta, = . The expec- 

tation value (■ ■ • )s is taken with respect to the distribu- 
tion of TiS. The conductance cumulants may be obtained 
from 



((.9")) ^ i-d)' 



.9>.(q) 



q=0 



(4) 



For large N it is natural to seek the CGF in the form of 
the expansion^ 



<^,(q)= ^i'^i^)^' 



(5) 



j = -oo 



In Ref.i2flit was found that in such an expansion, for the 
symmetry class A, ip^^ (q) is an odd (even) polynomial 
of degree 2 — j in q with j odd (even). For all the WD 

and BdG classes, in the region where ipi''\q) may be 
expanded in non- negative powers of L/£, it may be shown 
by induction from Eq. ljA2|l in Appendix A and from the 

definition of the CGF, Eq. ®, that (^i-'^(q) is of the form 
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We seek A^^' j, {L/£) in the form of a rational function 
of L/£. By using Eqs. ((SJ and 1A2|I it can be shown by 
induction that the further condition <^ s impUes 

that A\^^ (L/£) takes the form 

41...... (V^) = <^L,uSL/e)-' + o [{L/i)-^-^] . (7) 

Thus in the hmit N ^ oo, L/£ —> oo, s = constant we 
obtain for aU the BdG and WD classes the expansion 

1 2— J C30 n 
J — — oo n— 1 fci , . . — 1 i 

This is essentially an expansion in the inverse powers of 
I 



large bare conductance and it is expected to converge in 
the metallic regime 1/-/V s <C 1. We will show that, 
in the presence of the TR symmetry, a):_^ j, with n > 2 
actually vanish for the BdG classes CI and DHL Thus 
for the classes CI and Dili, the cumulants higher than 
the second contain no components that are analytic in s 
at s = 0. 



From Eq. (|A2|I in Appendix A one obtains the first of 

the coefficients a^^"* = — 1. For the other coefficients we 
obtain after a lengthy calculation an algebraic recursive 
equation (see Appendix A for more details). 
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Here we have 



a 



kp I J hp^ J ■ ■ ■ 1 ^Pri 



,0) 



for ni — 1, 
for Til = n. 



[-{-l ,kp^^^^-^ .... ,kp^ 

Moreover, we denote 

6I1 = 6l(-n - m), 6*2 = 6i(n - 2)61(2 - n - m), 
6I3 = 61(1 - n - m) 

with 



6l(n) 



1 for n> 0, 
otherwise 



Note also that X]/=o = '^ith fcp^ =1. For the higher 
cumulants this recursion relation is conveniently evalu- 



increasing n, and (3) for given m and n they may be 
calculated, e.g., in the order of increasing subindices. 
Equation Q implies 



i-drni y: 4T.\.,i^" 



It then follows that, instead of 

((5")) - s"+''-o-^-2^ 2<n< 1/s, 



(10) 



(11) 



that generalizes Eq. to the remaining six symmetry 
classes, we obtain the scaling 



„n+S„ 



0.2- 



2 <n < 1/s, 



(12) 



that is down by one power of s relative to the result of 
Altshuler, Kravtsov, and Lerner. 

Calculations based on a nonperturbative microscopic 
approach-^ suggest that there exist additional nonpertur- 



ated by a computer. Because of Eq. ®, the coefficients Native contributions that become significant for n > 1/s 



(m) 
^ki,...,k„ 



are invariant under the change of subindices. 
The coefficients are evaluated (1) in the order of decreas- 



and which lead to log-normal tails of the conductance 
distribution but are not describable within the one- 
ing m = 1,0,-1,..., (2) for a given m in the order of parameter scaling approach. Since we cannot produce 
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such nonperturbative components by our method we re- 
strict our considerations to conductance cumulants of the 
order n < 1/s. 

References and missed the cancehation of the 
terms of the order 0(s"~^) since the numerical prefac- 

tors corresponding to a^^ "\ were not evaluated. The 
cancellation of the leading contribution of Eq . ^ has al- 
ready been pointed out for n = 3 in Refs. iSCk .3 It and 
This cancellation is specific to quasi-one-dimensional 
geometryjAS For higher n the cancellation of the leading 
terms has been proved in Appendix B. 

I 



d-'{{9')) 



d-'m 



d-'m) 



For a normal-metal wire (m/ — 1), the term of the or- 
der O(s^) in the second cumulant and the first term of 
{{g^)) agree with Ref. Isil For the symmetry class A 
(mo = 2, 'nT±= 1), the subleading term of {{g^)) coincides 
with Ref. Hi 

For the BdG classes CI and Dili (mo = 2, m; = and 
mo = 2,m/ = 2) the coefficients of the form a^™-* , a^™]^,^ 
with m < — 1 equal zero. By a similar argument as in 
Appendix B one may show that this implies that all the 
factors of the form a^™'' with n > 3 and any m vanish. 

Thus in the case of the symmetry classes CI and Dili, 
the cumulants higher than the second do not contain 
components which are analytic in s at s = 0. For these 



B. Conductance cumulants 



For the numerical values of (((?")), the cases n — 1,2 
have been covered, e.g., in Ref. 31 for the WD classes 
and in Ref. El (n = 1) and Ref. |2J (n = 1,2) for the 
BdG classes. Our Eq. (O provides an efficient derivation 
of these results. As an application of Eq. ^ we give the 
results for n < 6, 



symmetry classes, the three lowest cumulants contain a 
contribution that is nonanalytic at s = 0.^^ We consider 
it likely that there exists a residual effect of this nonana- 
lytic behavior in the higher-order cumulants too so that 
these cumulants do not vanish, but we cannot prove this 
by our perturbative method. If such a residual effect ex- 
ists it is special to the cumulants higher than the second 
in the context of classes CI and Dili that the nonper- 
turbative components are actually the leading contribu- 
tions in the diffusive regime. This absence of an ana- 
lytical expression would indicate that the diagrammatic 
or semiclassical pictures'^^ do not lend themselves to the 
interpretation of these higher-order spectral correlations 



1 ^ (mo 2m;) ^ _ g^^^^^^ _ 2) ^ 4m,(m, - 2)]—^ 

s 3mo 45mQ 

-f {[-21mo + 6ml + 2m,(-5mo - ^mi + 16)](mo - 2) + 8m;(m; - l)(m; - 2)}-^^ + 0{s^), 

y45mn 



2s2 

■'0 



2 8(mo - 2)s ^. ^ , M 4s^ , 
' - [(15 + 19mo-69m0(mo-2) + 32?7i;(m; -2)]— ^ + 0{s^), 



Ibmo 315mg ' ' ' " 'M725mj^ 

8(mo-2)s2 ^ . 32s3 

= - \Lk 3 + 635 + 200 793mo - 751 117mO mo - 2 + 353 792mz m, - 2 — — — — ^ 
1485mg 638 512 875mQ 

+ {[-291 960 - 1 371 774mo + 554 337m^ + (2 303 952 - 1 429 012mo + 205 356m;)m;](mo - 2) 

+ 286 720mKm, - l)(m. - 2)} 837 625mg + 
512fmn — 2)s^ 

— J [(17 006 315 + 15 370 851mo - 62 563 249mO(mo - 2) + 29 630 464m/(m, - 2)] 



155 925m(S 



54 273 594 375m§ 



n 229 952fmn — 2"!^* 

^— ^ + [(55 989 824 345 + 45 233 187 091mo - 192 229 424 511mO(mo - 2) 

3 919 486 725mQ 

91 546 451 968m,(mz - 2)] , + 0{s^), 

510 443 155 096 875mn 



6 045 601 792(7710 - 2)5^ 
1 747 609 738 875m6 



6\ 





(13) 
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in the BdG wires with TR symmetry. 

The factor (mg— 2)/mo is a consequence of W(A)L. For 
the WD ensemble with broken TR symmetry (mo — 2, 
m; = 1) we obtain only even (odd) powers of s in ((fl")) 
with n even (odd), in agreement with Refs. js^ and Isil 
whereas for the symmetry classes AI, All, C, and D all 
the powers are present. 



The next correction terms are of the order ©(s^). Here 
we have adopted the notations 



Cl 



IqIq (?■ d {ma - 2mi)deVto 

— ' = ~r~^^ ^0 ^ 5 — I. ' 

e n s AniQn 

[(3too ^ 8mi){mQ — 2) + Ami{mi — 2)]deVt{)t 



(18) 



Current cumulants 



As a second application of Eq. we have calculated 
the current cumulants Cj familiar from the context of full 
counting statistics and defined^'^ in terms of the generat- 
ing function S{x) ^ - InEjv ^*o(^) exp{iNx)], 



(14) 



Here Ptg (iV) is the probability of N electrons traversing 
through the sample in a time to while x is a so-called 
counting field. Since we have 




E 



(m) - 

a) 's 



(15) 



Equations and (fT?^ provide a straightforward deriva- 
tion for the leading contributions of Cj since for these 
terms, in Eq. only the sums in square brackets 
with n^m = j — ni^pi — 1 contribute. For the lead- 
ing contributions of Cj, our results agree with those in 
Ref. 13 It can be shown by Eq. © that the correc- 
tion terms of certain order 0(s~™) depend similarly on 
Too and to; with all the j's (through co, ci, . . .) and dif- 
fer only by the numerical coefficients. This was proved 
in Ref. '35' for the terms of the order 0{s^). While the 
weak-localization corrections to the current cumulants 
have been studied before, e.g., in Ref. 36, to our knowl- 
edge, the exact numerical values for j > 2 have not been 
reported. Except for the first two cumulants, the ratio 
of the numerical factors before the universal correction 
term co and before the bare first cumulant Qo is for even 
cumulants of the order, but smaller than unity and for 
odd cumulants smaller by about a factor of ten. 



the equation^ 



III. DISCUSSION 



(16) 



yields for the current cumulants 



Cl =Qo + Co + Cl, 

^ 1 ^ 1 3 ^ 1 ^ 1 23 

Ct = — Un H Co Cl , Cq — (Jn H Cn H Ci , 

3^ 15 7 ' 15^ 315 105 ' 

n 1 n 11 401 

O4 = Un Co Cl , 

105^ 1575 3465 ' 

1 ^ 1 18 101 

C5=-— Qo-^co + :;,^ci, 



315 315 
47 221 24 433 

:Co 



:Ci, 



14 189 175 945 945 
27 ^ 811 62 993 



5005" 2 027 025 4 922 775 
_ _3_ 1 790 851 664 157 

^ ~ " 715*^° ^ 516 891 375*^° " 70 509 285*^^' 

233 98 299 813 1 095 799 



36 465 



206 239 658 625 
23 610 799 591 



969 969""" 3 781 060 408 125 
10 053 185 861 

^ 3 478 575 575 475^^' 



Co 



204 105 825 



(17) 



In conclusion, we have presented a recursion equa- 
tion covering the asymptotic behavior of the higher-order 
mesoscopic fluctuations for seven universality classes. We 
evaluated the values of the conductance cumulants (((?")) 
with n < 6 and the weak-localization corrections to the 
current cumulants Cj for j < 10. We discovered two 
qualitative results: (1) conductace cumulants of order 
larger than two and smaller than 1/s with s the length 
of the wire in units of the number of channels times mean 
free path scale with s with one less power than expected 
on the basis of naive scaling analysis and (2) the same 
cumulants are all vanishing in an expansion in powers of 
s in the two BdG symmetry classes characterized by TR 
symmetry, the fact from which we conjecture pure non- 
analytic dependence on s. As far as the noninteracting 
DMPK model is valid, P{g) deviates from the Gaussian 
shape only slightly. These deviations may, however, in 
principle, be detected by generating a large number of 
uncorrelated disorder realizations by repeatedly heating 
and cooling the samplei^ 
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APPENDIX A: DERIVATION OF THE 
RECURSION EQ. (9) 

The DMPK equation Q may be expressed through 
TiS instead of A^s (for the WD ensembles, see Ref. l29l) . 



Rewriting Eq. (Al) of Ref. |2J, which is a direct conse- 
quence of the DMPK equation, in terms of r^s yields for 
exp(— q • T) an evolution equation 



(moN + 1 + mi — Too) 
2N 



, OC r- y k~l \ y' k — 2 

^ fc=l L \ ^ ^ (=0 



klqkqi{Tk+i - Tk+i+i) + ^ kqk 



k,l=l 



fc=i 



1 - — I k(Tk+i - Tfe) + Tfe+i H Tk 



exp(-q • T) 



(mo,m;) 



(Al) 



For the CGF </j,(q), the terms of the order 0(iV"+i) 
imply 



J 



/II ^ OO ^ rfcl — 1 / 1 Q (j) c. (m— 7 + 1) 

mo („,.) 1 + to; -mp) , ^ jTOo / ' dipl 

^'^^^ ^ ^^^^ '^^^ £ ^''^l~^S vS^...-. -a^^ri^TT 

] mo ^ ^ _ (mo -2) ,^^ _ d^^\ ^ (TO, - 1) 9y<"'+^) 

2 ^ ^ 2 5g/c+i dqk J 2 Sg^i+i 



fc2 = l 

(too - 2to,) dip'. 



'fci+fc2+l 

(m+1) 



(A2) 



Making use of the expression for ips'{q) given by 
Eqs. © and JHl and considering the terms of the or- 
der 0(5^™ Jl"^^ gfc.) we arrive, after a lengthy algebra, 
at an algebraic recursive equation, our Eq. 

APPENDIX B: CANCELLATION OF THE 
TERMS 0{s"-'^) FOR {{g")) WITH 71 > 3 

Let us consider Eq. with a factor of the form 

Ogi^2,^---,9„ on the left-hand side (Ihs). On the right-hand 
side (rhs) of Eq. 0, only the sums in square brackets and 
the r sum on the third line may then contribute. There 
are no factors of the form CLqi^q2 ^^...,qj^ in Eq. ip. The fac- 
tors flgi^gj^ga vanish for all the WD classes 1^ Equation 
implies ag7,g2,«3 — fo'" the symmetry classes since the 
terms containing mi cannot contribute to agi,g2,«3- It is 



thus sufficient to show that on the rhs, the joint contri- 
bution of (i) sums containing products of factors of the 
form Uql'' and aqi^g"^...^<j„ , (ii) sums containing products 
of factors 0^2^53 and agi.(j"^...,g„_i , and (iii) factors of the 
form aqi^g"^...^(j„_i vanishes. 

For the induction argument, let us consider Eq. 
with TO = 2 — no, n = uq fixed and assume that we have 

flgi.92 °^-.9no-i — fo'' ^-ll positive integer values of g^. 
The induction assumption implies that the contributions 
of the forms (ii) and (iii) vanish such that we have 

ag?^2°^.,9„o = ^ "^^th qi ^ q2 = ■ ■ ■ ^ qno = I- This im- 
plies that also the contribution (i) vanishes and we have 

ag?,(j2 °^..<;nf, = with arbitrary qiS. The induction as- 
sumption holds with no = 4 since we have 0^1,^2,93 — 0- 
Thus we obtain a^qi,cll]...,q^ = for arbitrary n > 3. 
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